Microscopic understanding of the orbital splitting and its tuning at oxide interfaces 
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By means of a Wannier projection within the framework of density functional theory, we are 
able to identify the modified c-axis hopping and the energy mismatch between the cation bands 
as the main source of the ti g splitting around the F point for oxide heterostructures, excluding 
previously proposed mechanisms such as Jahn- Teller distortions or electric field asymmetries. Inter- 
facing LaAlOs, LaVOs, SrVOs and SrNbOs with SrTiOs we show how to tune this orbital splitting, 
designing heterostructures with more olxy electrons at the interface. Such an "orbital engineering" 
is the key for controlling the physical properties at the interface of oxide heterostructures. 

PACS numbers: 73.20.-r, 73.21.-b, 79.60.Jv 



When an atom is part of a periodic arrangement, such 
as a solid, the spherical symmetry of its potential is low- 
ered with respect to the case of a free atom. Typical is 
the example of transition metal ions surrounded by oxy- 
gen octahedra in cubic perovskites. Consequently, e.g., 
the five d orbitals split: In bulk SrTiOs one has three t2 g 
orbitals (d xy , d xz and d yz ) and two e g ones (d x i_ y i and 
d 3z 2_ r 2). The physics of transition metal oxides is deeply 
influenced by the further (finer) splittings of the t2 g or e g 
orbitals, for instance when distortions of the octahedra 
are energetically favored [l| . Both from the point of view 
of basic materials research and from that of technological 
development and device fabrication, it would be fascinat- 
ing if one could control such deviations from the perfect 
cubic perovskite structure by means of some external ad- 
justable parameters. Yet, this is something not easy to 
do in a flexible and controlled way in bulk materials. 

The recent breakthrough in growing oxide heterostruc- 
tures, such as LaAlOs grown on SrTiOs, offers a new 
possibility to tune the orbital degrees of freedom. The 
interface breaks the translational, and hence the cubic, 
symmetry. As a consequence, the three t2 g orbitals of 
the Ti atoms close to the interface split |2[. Evidences 
for similar effects come also from experiments on a bare 
SrTiOs (001) surface upon cleavage [!, 0]. The physical 
properties of the entire heterostructure, such as super- 
conductivity [5], phase separation [6], magnetism 7^9], 
etc., depend on the t2 g electrons at the interface [id. 
These effects have been revealed in a series of stunning 
experiments [l2|, [l3| , which also demonstrated that in- 
terfaces in layered oxide heterostructures can not only 
be seen as a way of tuning orbital splittings for specific 
desired purposes, but also for engendering new physical 
effects that are absent in the constituent bulk materials. 

A sine qua non condition for a successful "orbital engi- 
neering" in these layered systems is therefore the under- 
standing of the mechanism behind the above-mentioned 
lifting of the t2 g degeneracy. Several density func- 
tional theory (DFT) calculations have been applied to 



LaA10 3 /SrTi0 3 (LAO/STO) heterostructures [14H21]. 



They clearly show that the d xy band with Ti charac- 
ter close to the interface is no longer degenerate with 
the d xz and d yz bands. Some mechanisms for this ef- 
fect have been proposed such as Jahn- Teller distortions 
(l^ , a crystal field splitting because of the Sr 2+ vs. Ca 3+ 
asymmetry at the interface [17|, [l8| and a wedge-like or 
quantum- well- like potential at the interface [5[ [lo| . A 
simple, microscopic understanding of this mechanism is 
however hitherto lacking. 

In this Letter, we determine the relevant parameters 
governing the splitting between the d xy and d xz /d yz or- 
bitals at the interface by a Wannier function projection. 
We show that the c-axis hopping and cation energy mis- 
match is of primary importance, ruling out other poten- 
tial mechanisms. This understanding also allows for a 
tailor-made orbital splitting, which in turn controls the 
physical properties of oxide heterostructures. 

In the right panel of Fig. [I] we show the bands with t2 g 
character for bulk SrTiOs as calculated by Wien2k in the 
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FIG. 1: Right panel: Band structure (green lines) for bulk 
SrTiOs of two of the three t^ g bands: the d xy (red) and the 
dy Z (black). The thickness is proportional to the orbital char- 
acter. Left panel: Crystal structure of bulk SrTiOs and il- 
lustration of the maximally localized Wannier function with 
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paramagnetic phase employing the generalized gradient 



approximation (GGA) [22j. The bulk structure is per- 



fectly cubic, therefore the three bands are degenerate at 
the T-point. Along the T-X direction the xy and xz (not 
shown here) bands are strongly dispersive, due to the 
large hopping amplitude along the ^-direction, while the 
yz orbital (shown in the left part of the Figure) overlaps 
much less in that direction. The values of the hoppings 
are listed in the Table in the Supplementary Material. 

A very simple way of addressing the question of the ti Q 
splitting is to build a symmetric heterostructure with (for 
the case of LAO/STO) one AIO2 layer in the middle, two 
LaO layers around it and alternating TiC>2 and SrO layers 
on both sides (see Fig. [2j). Heterostructures of this kind 
have two n-type of interfaces: Ti02|LaO and LaO|Ti02- 
This implies that the interface will always be metallic, 
mimiking the situation of an n/p heterostructure above 
the critical thickness. The calculations have been done 



using the Wien2k code [22| with GGA and optimizing 
the internal coordinates. 

Several groups have analyzed the band structure of 
LAO/STO. The picture coming out of the DFT calcu- 
lation and supported by a number of experimental find- 
ings, suggests the presence of two types of carriers: a 
more itinerant one extending deep into the STO part 
and one more confined at the interface, primarily made 
of d xy electrons 0, Q . The former is easier to get trapped 
by impurities, defects or in the form of lattice polarons. 
Therefore, by tuning the splitting at the T point between 
the d xy and d yz bands we can control these two compo- 
nents and, in turn, the metallic or insulating behavior 
of different heterostructures. We consider here different 
heterostructures and for each of them we calculate the 
band structure focusing on the orbital character close to 
the interface. 

A very transparent description can be achieved by ex- 
tracting Wannier functions of d character and calculat- 
ing the tight-binding Hamiltonian arising from them. To 
our knowledge, this has been never done in the context 
of the heterostructures but it is actually a very power- 
ful and simple way of describing the electronic properties 
of the different cations at, as well as close to, the inter- 
face. In our case we consider maximally localized Wan- 
nier functions, because the resulting low-energy Hamil- 
tonian can then also be very naturally used to perform 
future many-body calculations. The Wannier projection 
was performed with the Wien2Wannier package [23| in- 
terfacing Wien2k to Wannier90 [24]. 

The two pieces of information that we are going to fo- 
cus on are i) the local term of the Hamiltonian, which 
gives the energy position of a given orbital (or, more pre- 
cisely of the center of gravity of the corresponding band) 
and ii) the nearest-neighbor hopping terms between dif- 
ferent t2 9 orbit als, parallel and perpendicular to the inter- 
face plane. From the projection onto maximally localized 
Wannier functions we directly get both. An analysis of 




FIG. 2: (a) and (b): Band structure for the LaA10 3 /SrTi0 3 
heterostructure with one AIO2 layer symmetrically sur- 
rounded by two LaO layers and STO. The thickness of the 
lines highlights the character of (a) the interface layer, (b) 
the second Ti02 layer. Shown are also the Wannier orbitals 
for the two layers, (c)-(e): Band structure for LaVOs/ (c), 
SrVOs/ (d) and SrNb0 3 /SrTi0 3 (e). The lowest parabolic 
bands (thin green line) are almost completely of V or Nb d xy 
character, panel (c)/(d) and (e), respectively. In all panels 
red (black) indicates xy (yz) character. 



the onsite energies for the different heterostructures al- 
lows us to identify the mechanism causing the shift of the 
different bands: the polar discontinuity between LaO^ + ) 
and SrO(°) layers or the presence of electrostatic poten- 
tial shift induced by the accumulation of charge at the 
interface, or possible octahedra distortions 0,l2l|. The 
hopping amplitudes between the Wannier orbitals give us 
information about the splitting at the T-point and about 
the degree of itinerancy of each orbital, another point of- 
ten discussed in these "subband" studies Q, Q . From our 
Wannier functions we can also analyze the localization of 
the orbitals in real-space as well as establish the amount 
of orbital reconstruction induced by the presence of the 
interface. 

In Fig. [2^l and [2)3 the band structure for the n/n 
LAO/STO heterostructure is shown. Even in the pres- 
ence of the interface, the Wannier functions are very close 
to ideal ti g orbitals. Therefore in Fig. [2] we only show 



the d yz one (the d xz is symmetrically related), for which 
a small degree of asymmetry can be observed. Unlike for 
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the bulk, we observe a splitting of the t^g bands at the in- 
terface, in agreement with earlier calculations [2|, |3j: The 
bottom of the d xy band orbital at the interface is lower 
than that of the lowest d yz . In addition to that, the bands 
with a component in the z direction (of which only the 
d yz is shown here) give rise to several branches close to 
the interface, which are due to the quantum confinement 
effect along the z direction. 

The major advantage of our Wannier projection is that 
we can now identify the mechanism behind the 250meV 



d yz splitting observed at the Y point. Indeed, 



from our single-particle Hamiltonian in the Wannier ba- 
sis, we can reliably extract the local terms s^ y and Sq Z . 
First, we can quantify the band bending at the interface: 
As one can see in Fig. [3b this amounts roughly to 0.3 
eV across three layers, in agreement with experimental 
indications [251 as well as with previous DFT calcula- 
Second, and even more important, since 
and Sq Z is as small as 50meV 




tions 

the difference between Sq V 



at the interface (see Supplementary Material) we can un- 
ambiguously rule out the polar discontinuity as the main 
source for the d xy vs. d yz splitting observed at the T 
point. 

A second effect determining the d xy 
at the T point is the following: When LAO is grown on 
STO a strong reduction of the vertical dy Z -dy Z and d xz - 
d xz hoppings occurs due to the presence of the insulating 
LAO overlayer (a similar mechanism has been already 
discussed for Ni-based heterostructures [I?}). The hop- 
ping amplitude along the z-axis involving the d xy orbital 
is instead small both for bulk STO and for the LAO/STO 
(see Table in the Supplementary Material). Since only 
the dyz and d xz hopping processes along z are affected 
the degeneracy with the d xy is lifted and the bottom of 
the dyz and d xz bands is shifted up (green/black band 
in Fig. [2ji). Quantitatively, this second effect coming 
from the hopping reduction dominates: While the dif- 
ference between the xy and the yz local terms is about 
50meV, the suppression of the hopping along the z direc- 
tion accounts for the remaining ~200meV inband split- 
ting, bringing the total splitting at the T point to the 
correct value which can be read off from Fig. [2Jl. 

Now that we have the two "control knobs" , the reduc- 
tion of c-axis hopping and the energy position of the t2 g 
bands of the different cations, we can use them to engi- 
neer the d xy vs. dyz splitting at the T point back to the 
vanishing bulk value or even to reverse it, pushing the d xy 
orbital above the d yz . Since the biggest effect is coming 
from the c-axis hopping reduction we first focus on the 
case of a heterostructure in which the polar discontinu- 
ity given by the overlayer is similar to that of LAO/STO 
but in which the hopping along the z direction is large. 
This way we expect to bring the bottom of the d yz band 
down in energy. Such a situation is achieved in the case 
of LaV0 3 /SrTi0 3 (LVO/STO), shown in Fig. [2fc. The 
bottom parabolic band (thin green line) is of V xy char- 
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FIG. 3: (Color online) Panel a: level scheme for SVO/STO for 
in-plane (V-V or Ti-Ti) and out-of-plane (V-Ti) hopping pro- 
cesses. Panel b and c: Layer dependence of the average over 
the three t<2 g bands of, respectively, the on-site energy and of 
the dy Z -dy Z hopping amplitude the along the z direction. The 
dyz splitting dashed lines represent the bulk values. 



acter. However, now the band with d yz character is a 
Ti-V hybrid. Like in LAO/STO, the LaO and V0 2 lay- 
ers in LVO/STO have, respectively, nominal charge +1 
and -1, but the presence of d electrons on V now lets the 
yz hopping survive. In Fig.[2k we indeed see that, at the 
interface, the d xy band has been pushed above the low- 
est (Ti-V hybridized) d yz branch at the T point, giving 
a situation totally different from LAO/STO. 

As next "orbital tweak" we want to push the d xy or- 
bital further up in energy in order to design an inter- 
face with dyz/d xz electrons only. This can be done by 
replacing La with Sr, i.e. by growing SrVOs/SrTi03 
(SVO/STO) [28|. Indeed, since the difference between 
the V and Ti Wannier local levels in the heterostucture 
is about 1.3eV (see Fig. [3)3 and Table in the Supple- 
mentary Material) we expect the situation illustrated in 
the scheme shown in Fig. [3^: The Ti d xy band is higher 
in energy and a mixed Ti/V band with d yz character is 
formed. Both effects can be observed in the actual calcu- 
lation shown in Fig.[2]i, where the Ti/V d yz band crossing 
the Fermi level is represented by the green and thicker 
black lines. Such Ti/V d yz band comes from the fact 
that, unlike the d xy orbital, the d yz and d xz orbitals of 
Ti and V hybridize to each other along the z-axis at the 
interface. The same hybridization is also responsible for 
the d xy vs. dyz splitting of the V states: The d yz orbital 
is much higher in Ti than in V and this "pushes" the V 
dyz orbital up in energy, as shown in the level scheme of 
Fig.Ek 

We have therefore designed a heterostructure with d yz 




Al/V/Nb Ti Ti Ti Ti 



Layer index 

FIG. 4: (Color online) Layer dependence of carrier densities 
of d xy (a) and d yz (b) orbitals 



(and d xz ) carriers only at the interface. To achieve this, 
the big energy mismatch between Ti and V t<i g orbitals in 
SVO/STO has been exploited. A case in which such mis- 
match is instead as small as 0.2eV is the SrNbOa/SrTiOa 
heterostructure (SNO/STO). As a result, the splitting 
around T looks close to that of an undistorted bulk per- 
ovskite with three degenerate t2 g bands. In other words 
with SNO/STO we can tune the t2 g splitting in such a 
wat that we get a seemingly "accidental" degeneracy as 
in the bulk case. Considering also the contribution from 
the hopping terms we can account for the ~0.7eV differ- 
ence in the position of the Nb d xy and Nb-Ti hybridized 
d yz bands observed in Fig. [2^. 

The ability of tuning the orbital splitting at the V point 
of oxide interfaces is not only per se intriguing but it also 
guides us to make heterostructures with desired prop- 
erties: insulating, metallic, with more or less localized 
carriers at the interface, etc. To illustrate this, we show 
the xy and yz occupancies in Fig. |4] for the four het- 
erostructures considered here (the behavior of the hop- 
ping terms and of the local levels is shown instead in 
Fig. [3b and c respectively and the values are reported 
in the Supplementary Materials). In the lower panel 
of Fig. E] the yz component is shown. This extends 
into the STO part and, since it has a large degree of 
itinerancy it is the component which gets more easily 
trapped hence it hardly contributes to transport (i.e. it 
is not observed in Hall measurements). The upper panel 
shows instead the xy electrons, a part of which forms 
a markedly two-dimentional electron liquid at the inter- 
face. Evidently, SVO/STO has no Ti xy carriers with 
high mobility. SNO/STO has a slightly larger fraction 
of d xy electrons and LAO/STO has more. On the ba- 
sis of this argument, the most metallic heterostructure 
should be LVO/STO. Yet, since V has two d-electrons 
in LaV03, strong correlation effects can contribute to 
localize them, thereby reducing the number of carriers 



available for transport. Experimentally, LAO/STO and 
LVO/STO with polar discontinuity have shown conduc- 
tivity with highly mobile carriers [HI, |2jJ [lo| . Our the- 
oretical study suggests that SNO/STO will have similar 
transport properties without polar discontuinity, whereas 
SVO/STO will have no xy carriers on the Ti interface 
layer. 

In conclusion, we have extracted maximally localized 
Wannier functions for four different heterostructures. 
This allows for a micoscopic understanding on how to 
tune the orbital degrees of freedom at the interface. The 
d xy vs. d yz splitting at the T point is mainly determined 
by the hopping amplitude along the z direction which 
is strongly influenced by the presence of the interface as 
well as by the energy difference between the Ti bands 
and those of the capping cations. We also showed that 
upon changing the material grown above STO we can 
tune this hopping and the shift of the local levels in such 
a way that the amount of d xy electrons can be engineered 
in a tailor-made way. This can be exploited to increase 
the amount of d xy electrons and, hence, to make more 
metallic heterostructures. 
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